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Abstract 

We describe how unintegrated parton distributions can be calculated from conventional 
integrated distributions. We extend and improve the 'last-step' evolution approach, and 
explain why doubly-unintegrated parton distributions are necessary. We generalise kf 
factorisation to (z, /ct)-factorisation. We apply the formalism to inclusive jet production 
in deep-inelastic scattering, mainly at leading order, but we also study the extension to 
next-to-leading order. We compare the predictions with recent HERA data. 



1 Introduction 



Conventionally, hard processes at proton colliders are described in terms of scale-dependent 
parton distributions a(x,/i 2 ), where a = xg or xq. These distributions correspond to the 
density of partons in the proton with (plus) momentum fraction 1 x, integrated over the parton 
transverse momentum up to k t = fi. They satisfy DGLAP evolution in the factorisation 
scale fi. However, for exclusive processes, parton distributions unintegrated over k t are more 
appropriate. The unintegrated distributions, f a (x, k 2 , /i 2 ), have the advantage that they exactly 
correspond to the quantity which enters the Feynman diagrams and therefore allow for the true 
kinematics of the process at small x, even at leading order (LO). We will explain how the exact 
kinematics may be restored for general values of x. 

-"■The plus and minus components of a parton with 4-momentum k are k^ = k°±k 3 . In the infinite momentum 
frame, the plus momentum fraction x = k + /p + becomes the longitudinal momentum fraction of a proton with 
4-momentum p. 



The distributions depend on two hard scales, k t and fi, and so the evolution is much more 
complicated. For example, the gluon distribution f g (x,k 2 , /i 2 ) satisfies the CCFM evolution 
equation pQ based on angular ordering of gluon emissions along the chain, in the approximation 
where only the 1/z and 1/(1 — z) singular terms of the splitting function P gg {z) are kept. So 
far, working with this equation has only proved possible with Monte Carlo generators [2]. 

However, in Ref. [SJ H] it was shown that it is possible to obtain the two-scale unintegrated 
distributions, f a (x, k 2 , /z 2 ), using single-scale evolution equations for h a (x,k 2 ) with the depen- 
dence on the second scale /i introduced only at the last step of the evolution. We call this 
the KMR procedure. 2 Two alternatives for the evolution of h a (x, k 2 ) were considered: (i) pure 
DGLAP evolution and (ii) a unified evolution equation |Hj which embodies both the leading 
log/c 2 (DGLAP) and logl/x (BFKL) effects, as well as including a major part of the sub- 
leading log 1/x contributions. As expected, the gluon and sea quark distributions, f a (x, k 2 , /x 2 ), 
extended into the k t > /i region, and indeed populated this domain more and more as x de- 
creased. An interesting result was that the unintegrated distributions obtained via the unified 
evolution of prescription (ii) were not very different from those based on the simpler DGLAP 
evolution of (i). It was concluded that the imposition of the angular-ordering constraint in the 
last step of the evolution was more important than including the BFKL effects. Here, we pay 
particular attention to probing the unintegrated quark distribution at larger values of x, so 
prescription (i) will certainly be a good approximation. 

We refine and extend the KMR last-step procedure [3] for determining the unintegrated 
parton distributions. First we note that in Ref. |3] angular ordering was imposed on both 
quark and gluon emissions; we correct this and only impose angular ordering on gluon emissions. 
Second, the KMR procedure was based on ^-factorisation [7j or the semihard approach ..$] (for 
a review, see Ref. (19 J in which the unintegrated parton distribution is convoluted with an 
off-shell partonic cross section where the incoming parton has virtuality — k 2 . This is only 
valid for gluons in the high-energy approximation where z — > 0, with z the fraction of the 
(plus) momentum of the parent parton carried by the unintegrated parton. Here, we generalise 
the notion of /^-factorisation and show that it is more accurate to calculate observables using 
doubly- unintegrated distributions f a (x, z,k 2 , fi 2 ), where the parton now has virtuality —k 2 /{I — 
z). 

In Section|H]we describe how the unintegrated parton distributions, f a (x, k 2 , /x 2 ), can be de- 
termined from the conventional integrated distributions a(x,/j, 2 ). Then in Section 0] we define 
the doubly-unintegrated distributions, f a (x, z, k\ , /i 2 ), and show how /^-factorisation is gener- 
alised to (z, /^-factorisation. The most direct way to test the unintegrated parton distributions 
is via inclusive jet production in deep- inelastic scattering (DIS). Inclusive jet production, par- 
ticularly in the current jet region, probes the unintegrated quark distribution in a similar way 
that inclusive DIS probes the integrated quark densities. The idea is that the LO diagram 
computed using (z, /c t )-factorisation will reproduce, to a good approximation, the results of 
the conventional LO QCD diagrams computed using collinear factorisation. This approximate 

2 An alternative formalism was given in Ref. 5 . 
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Inclusive jet production in DIS at LO 




Figure 1: A schematic diagram of inclusive jet production in DIS at LO which shows the 
approximate equality between, on the left-hand-side (a), the formalism based on the doubly- 
unintegrated quark distribution, and on the right-hand-side (b), the conventional QCD ap- 
proach using integrated parton densities, a(x,fi 2 ). 

equality is shown schematically in Fig. ^ The respective formalisms are presented in Section |S] 
and their predictions for inclusive jet production are compared with each other, and also with 
recent HERA data, in Section |H1 These sections not only compare the LO predictions, but also 
extend the comparisons to next-to- leading order (NLO). Section [7| contains our conclusions. 

2 Angular-ordered parton evolution 

We adopt a physical (axial) gauge, which sums over only the transverse gluon polarisations, so 
that the ladder-type diagrams dominate the evolution. Consider the evolution chain, illustrated 
in Fig. |2l simplified so that all the partons in the chain are gluons. The top of the diagram 
indicates some hard process with a factorisation scale \i. It is convenient to use a Sudakov 
decomposition of the 4-momenta of the propagator gluons: 

ki = Xip - faq' + k i± , (1) 

where p is the 4-momentum of the proton and q 1 = q + x B p is a light-like 4- vector. Here, q is 
the 4-momentum of the photon, with virtuality q 2 = —Q 2 , and x B = Q 2 /2p- q is the Bjorken-x 
variable. We neglect the proton mass since m 2 <Q 2 . In the Breit frame 

p={p + ,p-,p t ) = (Q/x B , 0,0), q' = {0,Q,O) and ^ = (0,0,^). (2) 

The emitted gluons along the chain have 4-momenta 

Pi = -h = fa-! -Xi)p+ (Pi - q' + p i± , (3) 

while the total 4-momentum going into the hard subprocess is 

P = k n + q = (x n - x B ) p + (1 - P n ) q' + k n± . (4) 
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P = k n +q 



Figure 2: Upper part of the evolution chain. 



Since the outgoing partons must be on-shell (p 2 = 0), we have 



(A - A 



p 



i-l, 



ar<_i(l - Zi) Q 2 ' 



(5) 



where Z; L = and the Sudakov (light-cone) variables of the propagator gluons obey the 

ordering: 

. . . > x n _i > x n > x B and . . . < /3 n _ x < /3 n < 1. (6) 

Colour coherence effects impose the angular ordering of the gluons emitted from the evolu- 
tion chain, originating from the destructive interference between the gluon emission amplitudes. 
The angle between the direction of the emitted gluons, with 4-momentum pi, and the proton 
beam direction must increase as we move towards the hard subprocess at the top of the evo- 
lution ladder. It is convenient to introduce a variable £j = p~ /pf. Then the rapidities of the 
emitted gluons are 

(7) 



Vi = - g lo §& = -log(tan(0i/2)), 



and the angular ordering, 



• • • < 0„-i <0 n <Q, 



(8) 



is equivalent to an ordering in & 



• • • < Cn-l < < 5, 



(9) 
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where 

S EE P-/P+ = il^L (10) 

provides the maximum allowed angle 6 via y/E = tan(9/2), assuming P 2 = 0. From (J3J) and 

c _ Pi _ f xbpv/q y = ( x B Pi y ^ 

where we have defined the rescaled transverse momenta 1 p i of the emitted gluons to be 

Pi = = — Qv^- (12) 

In angular-ordered (or CCFM) evolution, the factorisation scale \x plays the role of the maximum 
rescaled transverse momentum, so fi = x n Q\fE/x B . Therefore, the angular ordering (JSJ) can be 
written as 

. . . z n -ip n _ x < p n and z n p n < \i. (13) 



3 Unintegrated from integrated parton distributions 

It is informative to review how unintegrated parton distributions, f a (x,k 2 , /j, 2 ), may be calcu- 
lated from the conventional (integrated) parton densities, a(x,fi 2 ) = xg(x,fi 2 ) or xq(x,/i 2 ), 
in the case of pure DGLAP evolution. Recall that the number of partons in the proton with 
(plus) momentum fraction between x and x + dx and transverse momentum k t between zero 
and the factorisation scale /x is 

a(x,/i 2 ) — , (14) 

x 

whereas the number of partons with (plus) momentum fraction between x and x + dx and 
transverse momentum squared between k 2 and k 2 + dk 2 is 

f a (x,k 2 ^ 2 )^% (15) 

Thus the unintegrated distributions must satisfy the normalisation relation, 

rv 2 dp 

a(x,/i 2 )=/ — y f a (x, k 2 ,/2 2 ). (16) 
JO H 

The KMR proposal j3j to determine the unintegrated distributions was to relax the DGLAP 
strong ordering in the last evolution step only, that is, . . . k n -i,t <C k t ~ /i, where we have 
omitted the subscript n on the k t of the last propagator. This procedure is expected to account 
for the major part of the conventional next-to-leading logarithmic (NLL) terms, that is, terms 
like as{ots log/i 2 ) n_1 , compared to the usual leading logarithmic approximation (LLA) where 
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only terms like (as log /i 2 ) n are included. The procedure is as follows. We start from the LO 
DGLAP equation evaluated at a scale k t : 



da{x,k 2 ) = a s {k 2 ) 
d log k 2 2n 



E 

t>=g,q 



x 



dz Pah(z) b [ — , k t ) — a(x, k 



dccn a (0 



(17) 



where P a fo(^) are the unregulated LO DGLAP splitting kernels. The two terms on the right 
hand side correspond to real emission and virtual contributions respectively. The extra factor 
of ( in the virtual term avoids double- counting the s- and t-channel partons. The factor £ is 
equivalent to a factor of a half when integrating over ( and summing over b. 

The virtual (loop) contributions may be resummed to all orders by the Sudakov form factor, 



T a (k 2 , /i 2 ) = exp (- f^t ^ £ f <K C W 

\ Jk? K t Z7r l Jo 



(18) 



which gives the probability of evolving from a scale k t to a scale /i without parton emission. 
Differentiating, we obtain 



r, 



so that the DGLAP equation can be written in the form 



(19) 



da(x } k 2 ) = a s (k 2 ) 
d log k 2 2n 



E 



A- 2 



a(x,k 2 ) dT a {k 2 ,n 2 ^ 



T a (k 2 ,fi 2 ) dlogk 



(20) 



We define the unintegrated distribution to be 





fa(x,k t ,fl) = 



d log k 2 



[a(x,k 2 )T a (k 2 ^ 2 )} 



(21) 



This definition is meaningful for k t > /io, where /io ~ 1 GeV is the minimum scale for which 
DGLAP evolution of the conventional parton distributions, a(x, /i 2 ), is valid. Integrating over 
transverse momentum up to the factorisation scale we find that 



dA; 2 2 2 



Mo 



^(x,^ 2 )^^ 2 ,^ 2 )]^ 



= a(x,fi 2 ) - a(x,/io)T Q (//Q,/i 2 ). 
Thus, the normalisation condition (fTHj) will be exactly satisfied if we define 



72 /a^j , A* 2 ) 

kf 



\ a(x,/io)T a (/iQ,/i 2 ), 
Mo 



(22) 



(23) 
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so that the density of partons in the proton is constant for k t < //q at fixed x and fx. 

So far, we have ignored the singular behaviour of the unregularised splitting kernels, P a b(z), 
at z = 1, corresponding to soft gluon emission. These soft singularities cancel between the real 
and virtual parts of the DGLAP equation (fTTj) . After resumming the virtual part to all orders 
in the Sudakov factor (fTHj) the singularities must be regulated for the unintegrated distributions 
to be defined. The singularities indicate a physical effect that we have not yet accounted for. 
Here, it is the angular ordering caused by colour coherence, implying a cutoff on the splitting 
fraction z for those splitting kernels where a real gluon is emitted in the s-channel. 

We now apply the angular-ordering constraints of Section ^specifically to the last evolution 
step. For all other evolution steps, the strong ordering in transverse momentum automatically 
ensures angular ordering. The condition z n p n < fi (JTHJ) implies 

</x <{=> z<— V—, (24) 



1 - z fi + k t 

where, as before, we have dropped the subscript n specifying the last evolution step. Re- 
call, from the comment below (|T2*|) . that fi is entirely determined from the kinematics of the 
subprocess at the top of the evolution ladder: 



/i = Q-v/H = Q- MJL (25 ) 

X B X B W X i X B 1 

Equation (|2~l"j) applies only to those splitting functions in the real part of the DGLAP equation 
associated with gluon emission in the s-channel. By unitarity the same form of the cutoff must 
be chosen in the virtual part. We define £ max = 1 — C min = /•*/(/•* + ^t) and insert 0(Cmax — C) 
into the Sudakov factor for those splitting functions where a gluon is emitted in the s-channel 
and 0(C — Cmin) where a gluon is emitted in the t-channel. Note that there is no "coherence" 
effect for quark (fermion) emission and therefore the phase space available for quark emission 
is not restricted by the angular-ordering condition 



The precise expressions for the unintegrated quark and gluon distributions are 



f g (x, kl^) = T g (kl^) ^Sl J\z [p„(z) X -q (~ kt) (-t- - ^ + Pqg ( Z ) X -g (£, k\ 

"(26) 



and 



fi + k t 
(27) 



3 This is in contrast to Ref. 3 , where a cutoff on the splitting fraction was applied both to quark and gluon 
emissions. Also, in 0, the scale /x was treated as a free parameter, which was chosen to be the hard scale of 
the subprocess, or a combination of hard scales. Here we fix fi using 
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The exponent of the quark Sudakov factor can be simplified using the fact that P gq (l — C) 
P qq ((). Then 



dCC^«(C)+ / dCC^(C) = o 



Cm in 



Cmax p 1 

dCP qg (C)+ / dCP OT (C) 



so that 



Tq( k t,V ) = exp - 



^d/c? a s (Kp / ' Cma: 



fc? «? 2tt ./,, 



dC^«(C), 
(28) 



(29) 



Similarly, the exponent of the gluon Sudakov factor can be simplified by exploiting the sym- 
metry P qg (l — C) — Pqg{Q- We have 



£ jTdC C ^(0 = 2n F j\t l -P qg {Q = n F j\p qg (Q, 



(30) 



so that the gluon Sudakov factor is 



Tg(ktiV 2 ) = exp - 



" dn 2 t a s (Kl 



k? K t 



2tt 



d((P gg (()+n F J d(P q9 (C) 



(31) 



where n F is the active number of quark-antiquark flavours into which the gluon may split. 

It is important to note that the starting point of our derivation is the LO DGLAP equation 
(|17j) . with LO DGLAP splitting kernels and one-loop running coupling. Therefore, in order 
for the normalisation to be satisfied, it is essential that we use a LO parton set where 
the integrated parton distributions have been determined using the same splitting kernels and 
running coupling. In Ref . |3] , the MRST99 parton set [12] was used, which has been determined 
using NLO DGLAP splitting kernels and two-loop running coupling, therefore was found 
not to be satisfied. Also, in Ref. [3] the angular-ordering constraints were not correctly applied 
and the Sudakov factor T a (/ig,/i 2 ) was omitted from We have checked numerically that 

our refined prescription now gives the exact normalisation of (j!6|) . 



4 Calculating the cross section 

We have defined unintegrated parton distributions, f a (x, k 2 , /i 2 ), valid for all values of x for both 
the quark and gluon. This was done by assuming that the transverse momentum of the parton 
is generated entirely in the last evolution step and then imposing constraints from angular 
ordering to regulate the soft gluon singularities. It now remains to specify the prescription for 
calculating observables such as cross sections. 

The penultimate parton in the evolution chain has 4-momentum = xp/z. In the final 
evolution step, it splits into a parton with 4-momentum k n = k = x p — (3 q' + kj_ and an emitted 
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r 




Figure 3: Schematic picture of the k t - factorisation formula (J34|) 



parton of 4-momentum p n = k n _i — k n . The variable (3 is specified by the on-shell condition, 
p 2 n = 0, which gives 

Hence k 2 = —k 2 /(l — z). The rapidity of the emitted parton is 

rj*"* = - log & = - log (33) 
' 2 2 * x B z{3 1 ; 

In the small x regime, where gluons dominate, the main contribution comes from the z — > 
limit, where k ~ xp + A; 2 ~ — k\ and the emitted gluon has a large positive rapidity. In 
this case, observables can be calculated from the /^-factorisation prescription. For example, for 
deep-inelastic scattering, the cross sections for the scattering of a virtual photon with transverse 
or longitudinal polarisation can be expressed in the form 

41 = f~ r# /,(*,*?, M 2 ) (34) 

J X J Q K t 

see Fig. El At small x, we would expect that the leading log(l/a;) terms would need to be 
resummed. However, in Ref. [3 it was found that the unintegrated gluon based on a unified 
BFKL-DGLAP equation was very similar to the unintegrated gluon calculated purely from the 
DGLAP equation, as in Sectional 

In [H] the fc t -factorisation approach was used to calculate the unintegrated gluon contribution 
to the proton structure function F 2 (x B , Q 2 ). The unintegrated quark contribution was estimated 
in a rather ad hoc manner. In 0] the normal on-shell partonic cross section was evaluated with 
off-shell kinematics to estimate the cross section for prompt photon hadroproduction. Again, 
the ^-dependence of the hard-scattering coefficient was neglected. 
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4.1 Generalising /^-factorisation 



Clearly, it is desirable to formulate a more general prescription for the calculation of cross sec- 
tions using unintegrated parton distributions. This prescription should be valid for both quarks 
and gluons and without taking the limit z —>■ 0. The "partonic cross section" will necessar- 
ily have some ^-dependence, therefore we must consider parton distributions, f a (x, z, fef, /x 2 ), 
doubly-wamtegrated over both z and k 2 , satisfying the normalisation conditions 



J X 



dz f a (x, Z,k t ,fl) = f a (x, k t ,fl) 



and 



N 

Jx JO 



d/c 2 £ / __ _ j 2 2 



f a (x,z,k t ,(i ) = a(x,fi ). 



(35) 



(36) 



These normalisation conditions are only satisfied for fixed x and //, independent of the integra- 
tion variables z or k t . Apart from the angular-ordering constraints, the distributions may be 
obtained from ({21} : 



ft 1.2 2\ _ ^ /i2 2\ a s{kt ) 



2tt 



(37) 



The explicit forms, including the constraints, follow from ()26|) and (|27|): 

a 5 (/c 2 ) 



/ g (x, z,k t , [l ) — T q {k t , /X 



and 



2tt 



(38) 



/■ / 7 2 2\ _ /7 2 2\ a s{kt) 

Jg{X, Z, K t , fi ) — -Lg[K t , \1 ) 



2vr 



^p ot (,)^(^ 2 )+p m (,)^ (* *?) e( 



yU + fct 



(39) 



The universal factorisation formula involving these doubly-unintegrated distributions, anal- 
ogous to (jMJ), is 



a- 



a; 



dz 



1.2 
ft t 



fa(x,z,k 2 ,fi 2 ) cPp L (x, z, k 2 , fi' 2 ) 



(40) 



where £ are now the partonic cross sections for an incoming parton with (plus) momentum 
fraction x and transverse momentum k t , which has split from a parent parton with (plus) 
momentum fraction x/z and zero transverse momentum. We will refer to this generalised form 
of /^-factorisation as (z, fc^-factorisation. 

There will be an effective upper bound on the k t integration from kinematics, but note that 
there is no restriction to the domain kt < (i, as in conventional DGLAP calculations. For 
kt > /x, the Sudakov form factors T a (fc 2 ,/i 2 ) are defined to be 1. 
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Taking the limit z — > of crj|(x, 2, fc^, fi 2 ) in ()40|) we essentially recover the conventional k t - 
factorisation prescription of (|34j) . Alternatively, in the limit k t <C Q, we recover the conventional 
collinear factorisation prescription. 

Note that f a ( undefined for k t < fi ~ 1 GeV and also that (J23|) no longer 

applies since there is now a z-dependence involved. To approximate the kt < iio contribution 
of (J4"U|) . we choose to take the collinear limit k t <C Q in the hard-scattering coefficients, so that 



We then make the replacement 

'^dk 2 



k t <fJ.o K t <Q 



dz 



k 2 
ft t 



(41) 



(42) 



so that the (z, /ct)-factorisation formula (f^Oj) becomes 

f 1 f^dk 2 

a(x,i4)T a (iil,Li 2 )a^(x 7 Li 2 ) + / dz / — ± 



(7, 



(43) 

In the first term, the limit kt ^ Q must also be taken in the expressions determining a; and /j.. 
In the following, we will use (|4())l for brevity, with the understanding that the kt < liq region is 
to be dealt with as in (jjHj) - 



4.2 Motivation for the (z, /^-factorisation formula 

At this stage, it is perhaps unclear exactly how we should calculate the partonic cross sections, 
a^l(x, z, k 2 , fi 2 ), since the incoming parton is now off-shell with virtuality k 2 = —k 2 /(l — z), 
and so the usual /^-factorisation approach does not apply. This issue can be clarified by starting 
with the collinear factorisation formula one rung down. That is, 

b Jo \ x l z ) 

where we have chosen the factorisation scale to be k t , and b is the penultimate parton in the 
evolution chain of Fig. |21 so that a 7 b incorporates the last evolution step. From Fig. |21 we 
see that the parton b, with 4-momentum k n _i = xp/z, splits into a parton of type a with 
4- momentum k n = k = xp — (3 q' + k±, which then goes on to initiate the hard subprocess at a 
scale fi given by (J25|) . To derive formula (|4Uj) we need to show that the partonic cross section 
<T 7 * b can be factorised to give a partonic cross section for the 7*0 subprocess, cr 7 * a , with the 
remainder being absorbed into the definition of the doubly-unintegrated density, f a (x, z, k 2 , fi 2 ). 
This idea is illustrated in Fig. |U for the doubly-unintegrated quark distribution, and in Fig. |S] 
for the doubly-unintegrated gluon distribution. 
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Figure 4: Verification of (z, fct)-factorisation for the doubly-unintegrated quark distribution, 
f q (x, z,k 2 , fi 2 ), shown in the final diagram. In the first two diagrams the penultimate parton 
in the DGLAP evolution chain, with 4-momentum k n _x = xp/z, splits into a quark with 4- 
momentum k n = k = xp — [3 q' + k±. 




Figure 5: Verification of (z, fc t )-factorisation for the doubly-unintegrated gluon distribution, 
f g (x,z,k 2 ,fi 2 ), shown in the final diagram. In the first two diagrams the penultimate parton 
in the DGLAP evolution chain, with 4-momentum = xp/z, splits into a gluon with 4- 
momentum k n = k = x p — (3 q' + k±. 
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The squared matrix element can be factorised if we assume the LLA, so that only the 
leading 1/k 2 term is kept and terms not giving a logarithmic divergence in the collinear limit 
are neglected. We find that 

\Miy = 16^^^J2 P ^ \ M tV >< [1 + 0(0], (45) 

where \JvV a \ 2 represents the squared matrix element of the 7*0 subprocess, containing one 
power of as less than |.M 7 * b | 2 . We have used this method to derive the form of all four 
splitting kernels, P a b{z). It is crucial that we adopt a physical gauge for the gluon so that the 
splitting kernels are obtained from only the ladder-type diagrams. 

The extra terms of are proportional to (3 and so are negligible for either fc ( <C Q or 
z — > 0. Away from these limits, it is far from obvious that these "beyond LLA" terms will 
be small, a necessary condition for the factorisation to hold. We will observe that the main 
effect of the extra terms is to suppress the contribution from large z for gluon emission. In 
our approach, we achieve the same effect with angular ordering, so the extra terms may be 
neglected. 

\M 1 a \ 2 must also be evaluated in the LLA for the factorisation to hold, so terms of 
0(k 2 /Q 2 ) should be neglected when calculating this. This amounts to the replacement k — > xp 
in the numerator of \M. 1 a \ 2 , but not in the propagator virtualities in the denominator. Of 
course, x may have some kt dependence from kinematics, so some terms beyond the LLA are 
included in this respect. 

The phase space d<£> 7 * fe can be factorised easily to give the phase space d$ 7 * a : 

d 3 p n 



d$ 7 * b = d$ 7 * a 



2pl (2tt) 
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d$7 ' a (2^ d4 ^ ( ^ } (46) 
d$ 7 * a -—dxdk? 



167T 2 * x(l — z) ' 

where we have used d 4 A; = p ■ q dx d/3 d 2 k t and d 2 k t = k t dk t dip = Tfdk 2 , after integrating over 
the azimuthal angle <fi. The (3 integration absorbs the delta function, determining (3 as given 
by fl22). 

The partonic flux factor F 7 * a is not well defined since the parton a is off-shell and non- 
collinear with the photon. As in conventional fc r factorisation, we define it to be 4 

F 7*a = z p^b = zA . g = Axp . g ( 4? ) 

Finally, we have the relationship 

X Kt ATI * — ' 



4 Choosing another definition for the flux factor is a NLL effect. 
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To calculate the hadronic cross section, we insert ()48|) into (jUj) 

6 * a 

a 1 

where in the last step we recognise the "real" part of the doubly-unintegrated distribution given 
in (|37p. The (z, /^-factorisation formula ()40|) follows easily. 



5 Application to inclusive jet production in DIS 

The simplest process that we can consider to illustrate the use of the doubly-unintegrated 
parton distributions is current jet production in DIS. The subprocess is simply 7*g — > q at the 
top of the evolution chain. In the normal collinear factorisation approach, this diagram gives 
the parton model prediction for the structure function F 2 (x B ,Q 2 ). Indeed, measurements of 
F 2 (x B , Q 2 ) are used to determine the integrated quark distribution q(x, fi 2 ). In the new (z, k t )- 
factorisation framework of Section 0J where the incoming quark has transverse momentum k t , 
we produce a current jet with transverse momentum fc t and transverse energy Et = kt- The 
parton emitted in the last evolution step will emerge with transverse momentum — k t and 
transverse energy Et = k t . 

The inclusive jet cross section counts all jets passing the required cuts. The cross section, 
integrated over bins in y, Q 2 , E T and r\ is 

'Get) = / dy dQ 2 — ™ (1 + (1 -2/) 2 K P + 2(1 -yW p \ (50) 

X ^ ( E T - ^Tmin) © (#T ma x ~ E T ) 9 (rj - T] min ) 6 (r? max - T]) , 
jets 

where y = Q 2 / (x B s), and where the sum is over all jets with transverse energy Et and rapidity 
r\. The differential cross sections are easily obtained by dividing by the size of the bin, for 
example, 

: ff (jet)/(£Tmax - ^Tmin)- (51) 



dEj 

In Section |3] we gave the general prescription for calculating the cross section. Recall that it 
was necessary to consider the doubly-unintegrated parton distributions, f a (x, z, k 2 , /i 2 ), to keep 
the precise kinematics in the subprocess, without taking the limit z — > 0. We now check that 
this prescription reproduces with good accuracy the conventional LO QCD calculation with 
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integrated partons, where all 0(as) diagrams are included, not just the ones which give the 
leading dk 2 /^ term. With the (z, fc t )-factorisation approach, in addition to the jets produced 
in the hard subprocess, we must also count the parton emitted in the last evolution step with 
transverse energy Et = k t and rapidity given by (|33p. 

We also explain how the prescription may be extended to higher orders in perturbation 
theory. The conventional NLO QCD diagrams are at 0(a 2 s ). These include all real and virtual 
0{a,s) corrections to the LO QCD diagrams. The hard-scattering coefficients obtained from 
these diagrams are convoluted with NLO integrated partons, a(x,fi 2 ), satisfying the DGLAP 
equation with two-loop as and splitting kernels. Several codes are available which include these 
NLO QCD calculations. There is no longer a one-to-one correspondence between partons and 
jets. The 4-momenta of the outgoing partons should be passed through a jet algorithm to assign 
the partons to jets. At NLO in the (z, fc t )-factorisation approach, we continue to use the LO 
doubly-unintegrated partons constructed in Section El and only calculate the 0(acs) diagrams 
expected to dominate. 



5.1 Collinear factorisation approach at LO 

In the collinear approximation, the LO QCD Feynman diagrams are at 0{ots)- These are the 
boson-gluon fusion process, 7*g — > qq, and the QCD Compton process, 7*g — > qg, illustrated 
in Fig. HJb). These partonic processes give rise to two jets with equal transverse energy and 
opposite transverse momentum. There is a one-to-one correspondence between partons and 
jets. There are no singularities to be regulated and no cutoff is imposed on gluon emission. 

We now explain a few of the details involved since this calculation offers valuable insights 
into the (z, kt) -factorisation approach. The cut diagrams are illustrated in Fig. |H1 Note that 
the direction of fermion number flow is not indicated in these diagrams. The arrows indicate 
only the direction of the labelled 4-momentum and this is taken to be the same for both 
quarks and antiquarks. The contribution from diagrams (a) to (f) to cr rp need to be added 
together. Diagrams (a) to (d) have the same kinematics, so we calculate them first. We label 
the 4-momenta by 

x 

q = q' — x B p, l = —p, k = x p — (3 q' + k±, (52) 

x 

ji = k + q = (x - x B )p + (1 - (3) q' + k± and j 2 = I - k = -(1 - z) p + j3 q' - k±, 

z 

with x > x B . The 2-body phase space is 

i ^ = ( 2,r^Hi + ,-j 1 -h)-§^s(^^Ui) = 0%VO 2 2 ). (53) 

The two delta functions can be used to determine (3 and x: 

f3= x B zr aM = x B ( l _ z + r± ^ {l _ z + r y_ 4rz(1 _ A (54) 
x 1 — z) 2(1 — z) \ / 
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Figure 6: Cut diagrams contributing to inclusive jet production in LO QCD. 
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where r = k 2 /Q 2 . The flux factor is F 1 * a = 41 ■ q, so that 

_ dkj fxA 2 y, Z 2 1 

F~f* a 16tt \Q 2 ) x 2 {l-z)l-x B f3/x' [ } 

In practice, the condition x > x B ensures that only the x = x + solution contributes. 
The squared matrix elements of all six diagrams can be written in the form 

\Mf£\ 2 = \e 2 g 2 W» e,(q, X)el(q, A), (56) 

where A is either T or L and the initial factor of 1/2 is to average over the helicity of the incoming 
parton. Appropriate scales have been chosen for the two running couplings, e 2 = 4ira(Q 2 ) and 
g 2 = Anas{k 2 ). We have 

(a) M** = (Y,el)T R ^TT{ltYhrhliY]d P Al), (57) 



(b) 



q 

1 



(c) MT = e 2 C F -jT[^Yfrh u hY}d P Aj2), (59) 

(d) M*" = elC F ^-^^Yh V {Hi)l^il"]dAh), (60) 

where the colour factors are T R — 1/2 and Cp = 4/3. In an axial gluon gauge with a light-like 
gauge fixing vector n — q', 

d pa i!) = -g pa + n ^ + l ^ . (61) 
For diagrams (e) and (f ) of Fig. El the 4- momenta can be parameterised as 

q = q'-x B p, l = Xp, ji = £p + bq' + k± (62) 
and j 2 = I + q - ji = (X - x B - £) p + (1 - b) q' - kj_. 

with < £ < X — x B < 1 and < b < 1. This time the 2-body phase space determines 

b = x B rji and £± = l - \x - x B ± y/(X - x B )(X - (1 + 4r)x s )} . (63) 
Dividing the phase space by the flux factor gives 

a 2 ) ^ xl \i-ux-x R )iE\ ' (64) 



d^* q dk 2 fx, 



F-r*o 16tt \Q 2 J ££XZ\l-b(X-x B )/£\ 
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and the squared matrix elements are 

(e) M** = e 2 q C F + * ]_ Tc \tf{j x - 4)YJiY(f+ 4)Y] d pa (h), (65) 

(0 = e\C F {{ + * * Tr WQ/+ ki\H M d^h). (66) 

Averaging over the transverse photon polarisations in (|56j) . we have 



e M (g,T)e:(g,T)^-^, (67) 



while demanding that the longitudinal polarisation vector is normalised, [e(q, L)] 2 = 1, and 
satisfies the Lorentz condition, q ■ e(q, L) = 0, leads to 

^(Q, l ) = 7^( 2x ^ + ( 68 ) 

Gauge invariance ensures that the g M term does not contribute to the squared matrix element 
if all diagrams are included, courtesy of the Ward identity: 

q^M^ = = q u M» u . (69) 

Therefore, we are free to neglect the q^ term of (}68|) from the outset, so that 



Ax 2 

e^q,L)el{q,L)^-^ PllPv . (70) 



Finally, the contribution to the 7*p cross section from Fig. El (a), (b), (c) and (d) is 

4vr 2 ^ I" 1 [°°dk 2 ^ x B /x as(g) . . 

* T * = Z-V-J. dZ Jo 2. (71) 

x {P ?9 (z) -<?(-, // 2 ) fe + C^J ~?(~^ 2 ) fe + C^]}, 

while the contribution from diagrams (e) and (f) is 

A ^ ae l f\ v f°°dk 2 x B k 2 t ^ x B H a s (k 2 t 



1 



*v_ST ^ < (ay r dk * XB k ' X °/t as ^ k t) r Yn(Y ii 2 



(72) 
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where the coefficients are 

_ g {x + 2x B z-Axz-2x B z 2 + Axz 2 ) _ 8 (3 (1 - (3) x B (1 - z) z 

~ x (l-2z + 2z 2 ) ' ~ x (1 -2z + 2z 2 ) ' 

C T = A(x-2(3x-2/3x B z-2xz + A(3xz), C\ = 8 A(3 (1 - /?) x B z , 

c _ /5 (x + i B 2-iz-s i! 2 2 + 2i2 2 ) _Af3(l-f3)x B (l-z)z 

X (1 + Z 2 ) ' ~ X (1 + z 2 ) ' 

, _ -/3xz(l-z) , _ 

Ct " (x- Xb z) (i + z 2 y Cl -°> (73) 



(X-x B ) (X-x B -0 (^+(1-6) x fl )" 



where 



A = (5{l-z) [(x + x B z- I3x b z-xz){1 -2z + 2z 2 )] . (74) 



Note that for high Et jet production in LO QCD there are no infrared singularities from 
either on-shell propagators or soft gluon emission. We will take the factorisation scale to be 
fi = Et = kt, in order to compare directly with the approach based on unintegrated partons. 
The inclusive jet cross section calculated using (fTTj) and (|72j) was found to be in excellent 
agreement with the LO QCD predictions of the JetViP fi] and DISENT ^2] programs. 

At this point it is an interesting check to take the DGLAP limit, so that we insert 0(/z — k t ) 
and take the limit k t <C Q, so that the only contributions come from the ladder- type diagrams 
of Fig. El (a) and (c), and 



op = 



4vr 2 ae 2 f 1 f^dk? a s (k?) r_ , .x ,x 2 . _, , , x ,x 2x l 



1 

with x = x B and al p = 0. At lowest order 

Q 2 



F 2 (x B} ^) = ^(ap + ap)=J2e 2 q xq(x^ 2 ), (76) 



q 



leading to the well-known logarithmic scaling violation of F 2 , or equivalently the "real" part of 
the DGLAP equation for the (integrated) quark distribution: 

where the conventional choice of scale is /i = Q. To obtain the "virtual" part of the DGLAP 
equation for the quark distribution, it is necessary to consider loop corrections to 7*g — > q. Of 
course, for high Et jet production, it is not appropriate to take the limit k t <C Q. 
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Let us anticipate how this calculation would be treated in terms of unintegrated partons, 
where we would want to factor out the emission with 4-momentum j 2 in Fig. El (a) and (c) into 
the doubly-unintegrated quark distribution, f q (x,z,k 2 ,fi 2 ). For this to be possible, we must 
assume that Cj. = 1 = C£, and neglect all other contributions. The diagrams in Fig. El (d), (e) 
and (f) come from the subprocess 7*g — > qg, where the gluon is radiated off the final quark 
line. Such diagrams are strongly suppressed in an axial gluon gauge, due to one or more of 
the propagators having very large virtualities, and can be neglected. Similarly, for the crossed 
quark box diagram of Fig. El (b). Numerically, the terms proportional to (3 in diagrams (a) and 
(c) are found to be very small. The one exception is the term proportional to f3 in Cj>. This 
is negative and increasingly important as z increases, that is, it is a destructive interference 
term. In the case of our doubly-unintegrated quark, the same effect is obtained with an explicit 
constraint from angular ordering, so the term proportional to (3 is redundant. 

Ultimately, we will need to resort to explicit numerical comparison of (z, /^-factorisation 
with the conventional collinear factorisation approach in order to demonstrate the approximate 
equivalence of the two methods. 



5.2 (z, /ct)-factorisation approach at LO 

With the new (z, fcf)-factorisation framework developed in Section 0] the LO diagram is simply 
7*g — > q, illustrated in Fig.^a), where the incoming quark has 4-momentum k = xp — j3 q' + k±. 
The partonic cross section contained in (J4(Jj) is 

da£«(x,z,kl» 2 ) = d$^\Miy / F^«, (78) 
where F 1 q = Axp ■ q = 2xQ 2 /x B . Labelling the current jet by 

P = k + q = (x-x B )p + {l- fi)q' + k ± , (79) 
where x > x B , the 1-body phase space is 

= (2tt) 4 5 (4) (k + q - P) 5(P 2 ) = 2n5(P 2 ) 

{2ix) A 

where x± is given by with r = k 2 /Q 2 . Again, the condition x > x B means that only the 
x = x + solution contributes. The rapidity of the current jet in the Breit frame is 

V f eit = - log — = - log X/Xfl ~ 1 . (81) 
ip 2 B P-2 B l-f3 K J 

The squared matrix element, given by the cut diagram of Fig. [TJ is 

l-MS| 2 = ^Tr[^ 7 ^(^ + e,(q,X)e:(q,X), (82) 
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Figure 7: Cut diagram contributing at LO in the (z, fc t )-factorisation approach. 



where A is either T or L. We use the same formulae, (JoTj) and (J7UJ), to sum over the photon 
polarisations as before. 

Note that our approach is not gauge invariant since we do not include the complete set 
of cut diagrams shown in Fig. El Rather, we only keep the leading dk 2 / k 2 term coming from 
Fig. El ( a ) and (c). We rely on using a physical gluon gauge where the neglected diagrams are 
suppressed. We represent this approach by Fig. where the incoming quark is off-shell with 
virtuality —k 2 /(l — z). Strictly speaking, the Ward identity (JoTJJ) does not apply to Fig. [7| 
For example, the term of the longitudinal photon polarisation vector ()68|) gives rise to large 
cancellations between the contributions from Fig. El (a) and (b) to ensure that the Ward identity 
is satisfied. When the diagram of Fig. Efb) is neglected, as in Fig. the g M term in e^(q,L) 
gives a much too large o^. Therefore, we should not include the q^ term in e^(g,L); this is 
equivalent to an appropriate choice for the photon gauge. 

According to the prescription given in Section 0] we should only keep the leading dk 2 / k 2 
term in the squared matrix element and so terms explicitly of 0(k 2 /Q 2 ) should be neglected 
when calculating |A^|| 2 . This amounts to the substitution k = xp in the trace (jH^jl . leading 
to 

\Ml* q \ 2 = 4nae 2 g Q 2 — and \M£ q \ 2 = 0. (83) 

X B 

The partonic cross sections are then 

al q (x,z,k 2 ,fx 2 ) = q 2 - _ x B ^ x $(x-x+)e 2 g and a\ q (x, z, k 2 , /i 2 ) = 0. (84) 

Inserting into (J4(Jj) we obtain the hadronic cross section 

y * v An 2 a f 1 . f°°dk 2 x B lx v-^ 2 * < 12 2\ ,„ r < 
Q Jx Jo M l-x B f3/x ^ 

with x = x + . Again, it is an interesting check to take the collinear limit, k t <C Q, so that we 
insert 0(/z — k t ) and take \x = Q. Then, x — > x B , (3 — > and by the normalisation condition 
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we recover the parton model prediction for the proton structure function F 2 = Ft + Fl' 



F 2 (x B ,Q 2 ) 



Q 2 



(a T * p + a£ p ) = ^e 2 q x B q{x B ,Q 2 ). 



(86) 



47r 2 a 



Alternatively, taking the limit z — » of x and f3 in ()85|). then using the normalisation (|35|). 
gives a ^-factorisation prediction: 



with a; = x B (l + k 2 /Q 2 ). 

To test the assertion that the angular-ordering constraint mimics the major neglected terms 
in the LO QCD calculation of Section UTTl we can replace P qg (z) by P qg (z) {C a + C b ) and P qq (z) 
by Pqq(z) C c in the real part of the doubly- unintegrated quark, where the coefficients were given 
in (J73j) . The inclusive jet cross section calculated in this manner, with separate coefficients 
for the T and L contributions, is found to be almost unchanged, providing evidence that the 
destructive interference terms in the conventional LO QCD calculation have much the same 
effect as an explicit angular-ordering constraint. 

5.3 Towards a NLO (z, /^-factorisation approach 

It is beyond the scope of this work to perform a full NLO calculation within the framework 
of (z, /^-factorisation. Rather, at this exploratory stage, we aim to produce a simplified de- 
scription using the LO doubly-unintegrated partons and computing only the 0(ots) diagrams 
expected to be dominant. We do not want to include diagrams involving divergences which 
cannot be regulated by putting a cutoff on soft gluon emission from angular ordering. The ma- 
jor loop corrections are already accounted for by the Sudakov form factor (|18|). The diagram 
where a gluon is radiated from the final quark line is strongly suppressed in a physical gauge. 
This leaves the cut diagrams of Fig. |H]as the only contributions which should be included. It is 
debatable whether or not the crossed box diagram of Fig. |B(b) should be included. We choose 
to include it, although it gives only a relatively small contribution to the cross section. 

All diagrams in Fig. |H] have the same kinematics. An initial parton, with 4-momentum 
k = xp — (3 q' + k±, splits to a quark with 4-momentum k' = x' p — (3' q' + k' ± , which goes on to 
interact with the photon. The outgoing partons have 4-momentum j\ = k' + q and j% = k — k' 
where x B < x' < x < 1 and < (5 < f3' < 1. Including the parton emitted in the last evolution 
step, the Kt Jet package [IS] was used to cluster the three outgoing partons into jets using the 
inclusive k± algorithm in the Breit frame. Note that the diagrams of Fig. |S] naturally include 
the LO contribution of Fig. [7|in the limit that k t <C k' t . Therefore, the LO contribution does 
not have to be added in explicitly. 




(87) 
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We find that the 2-body phase space divided by the flux factor is given by 

( J) 2 E \ x - x ^ - 0- - ^ x> - (* - x *)P\~ x > ( 88 ) 



FT* a ~~ 167T.T 



where (3' = (3 + (x B R)/(x — x') and 
1 

2(W) 



x' ± = 1 m <!x(l -/?) +x B {l -(3-R) +x B r' (89) 



±y Ml -(3 + R)-x{l- (3)] 2 + x B r' [x B r' - 2 (z(l - (3) - x B (l -(3- R))\ 

with r' = k' 2 /Q 2 and R = \k t - k' t \ 2 /Q 2 . 

The cut diagrams representing the squared matrix elements are shown in Fig. |H1 Again, we 
write 

\Ml* a L \ 2 = -e 2 g 2 M*» e^q, A)<(ff, A), (90) 
where A is either T or L and the initial factor of 1/2 is to average over the helicity of the 
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incoming parton. We take e 2 = 4ira(Q 2 ) and g 2 = AnasifJ^), with /i^ = max(k t , k' t ). We have 
(a) = (j24)T R ^TTWYhl^'rjiY]d P M, (91) 



(b) W» = (J2 {k > + 1 - ky Tl W^^W + * ~ WW] d P *(k), (92) 



(c) = e 2 q C F ^iWYH¥l u hl"]d P .{j2). (93) 

In order to keep only the leading dk 2 /k 2 term, we make the replacement k —>■ xp in the 
numerator of these expressions, but not in the virtualities in the denominator. Inserting the 
partonic cross sections into (J4(Jj) we finally obtain 



q 

x 

x'=x 

X 



■ -" , 

Q \Q J Jo Jx Jo kt Jo X 27T 



|x-x fl /3- (1 -P)xf- {x-x B )P'\ 1 (94) 



{T fi / 9 (x, z, k 2 , /i 2 ) [C£ £ + 4 >£ ] + C F f q (x, z, k 2 , /i 2 ) C C T>L ) , 

where the coefficients are 

(1 - 2 /?' (1 - /?')) & (z' - ar B ) + (/?' (x fl - 2 ar') + x') ((1 - /?') x B + (2 /?' - 1) x') 



x ((!-/?') x s -x') 2 



4 (1-/5') p'x B {x 1 -x B +f3' {x + x B -2x')) 



a ((!-/?') x B -x') 



2 



b (i _ g) (^-^((1-2^^2 (g (2 x' - x B ) - xQ) 



r 6 



X 1 


- (i - 


/?') x fl ) 


((1 + /3-/3') (x-< 


) + (i + p-P'+r) x B y 








8(1 -/?') 2 /?'x s (x'~ 


x B ) 


X ( 


[x> - (1 - 


/?') x fl ) 


{(l + P-F) {x-x') 


1 + (1 + p - P' + R) x B ) ' 


(1 


-2/5' (1 


- /?')) a: 


(x'-x s )+x' ((2/3' 


-1) x fl + (1 - 2/3 /2 ) x') 



(x-x') ((1-/3') x B -x') 2 



4 (1-/3') P' 2 x B 



L ((1-/3') x B -x') 2 ' 

Inspection of the coefficient C T reveals a pole at z' = x'/x = 1, corresponding to soft gluon 
emission. We can regulate this singularity by appealing to angular ordering. The rapidity 
of the gluon, with 4-momentum j 2 , should be greater than the rapidity of the quark, with 
4- momentum ji. 



7 Brei t> Breit ^ /< ^ , = Q *\ J_P 

J J H' + \k t -k' t \ x B W x'/Xs - 1 
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Figure 9: Feynman diagrams contributing at "NNLO" in the (z, kt) -factorisation approach. 



This condition applies only to the diagram where a quark radiates a gluon, Fig. E^c), but not 
to the diagrams where a gluon radiates a quark, Fig. El (a) and (b). 



5.4 An estimate of the NNLO contribution 

The next-to-next-to-leading order diagrams have not yet been calculated in the collinear ap- 
proximation (NNLO QCD). As we explain later, the "NLO" calculation of Section 15.31 gives 
reasonable agreement with conventional NLO QCD. It is possible that a simplified "NNLO" 
(z, kt )-factorisation calculation may provide an estimate of whether the NNLO QCD corrections 
are likely to be important, especially at low Et and low Q 2 in the forward region, where there 
is a discrepancy between NLO QCD and the data. 

The four contributing diagrams, all of which have the same kinematics (phase space), are 
shown in Fig. El Diagrams (a) and (b) are the doubly-unintegrated quark contribution, while 
diagrams (c) and (d) are the doubly-unintegrated gluon contribution. Encouraged by the fact 
that the crossed quark box of Fig. Efb) gave only a small contribution, we may neglect the 
interference cut graphs arising from Fig. El as a first approximation, leaving only four squared 
matrix elements to be calculated. 

Our simplified approach provides an approximation of QCD, in which only ladder-type 
diagrams remain. The soft gluon singularities are regulated by angular ordering. There are no 
infrared singularities remaining. We can add an arbitrary number of rungs to the ladder and 
the answer will be finite. However, with more rungs, the number of neglected interference terms 
grows; it is likely that the approximate treatment of these terms by imposing angular-ordering 
constraints will spoil the accuracy of the method if too many rungs are added. 
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6 Description of HERA inclusive jet production data 



HERA data are available for inclusive jet production in DIS. We may therefore check how 
well the simpler (z, fc t )-factorisation approach is able to reproduce the conventional collinear 
factorisation approach, and at the same time see how well these calculations describe the data. 

Recall from Section El that at LO the (z, ^-factorisation approach is based on the simple 
7**7 ~~ > 1 subprocess driven by the doubly-unintegrated quark distribution, f q (x, z, k 2 , /u, 2 ), re- 
taining the full kinematics. On the other hand, in the LO QCD description the subprocesses 
are 7*g — > qq and 7*g — > gq evaluated with collinear kinematics and conventional integrated 
distributions, g(x,Q 2 ) and q(x,Q 2 ). 



6.1 Comparison with ZEUS data at high Q 2 

We now compare our predictions to the experimental data obtained by the ZEUS Collaboration 
[T3] . This data was taken during 1996 and 1997, when HERA collided protons of energy 
E p = 820 GeV with positrons of energy E e = 27.5 GeV at a centre-of-mass energy of y/s = 
^jAE p E e ~ 300 GeV. Rather than make cuts on the variable y = Q 2 / (x B s), ZEUS make cuts on 
cos 7, one of the angles used in reconstructing the kinematical variables using the double-angle 
method, where 

x B (l-y)E p -yE e 

cos 7 — — i ttt — • (97) 

x B (l-y)E p + yE e v ; 

In the parton model, the angle 7 corresponds to the direction of the scattered quark. In (|5T)|) 
we therefore set y min = and 7/ max = 1 and demand instead that cos 7 satisfies the ZEUS 
experimental cuts, —0.7 < cos 7 < 0.5. 

In Fig. El we show the rapidity distribution, dcr/d?7 Breit , integrated over Q 2 from 125 to 
10 5 GeV 2 and over Et from 8 to 100 GeV. The parton-to-hadron correction factors given in 
Table 3 of the ZEUS paper [T3] have been applied to the theory predictions. For the results 
presented, we used the MRST2001 LO parton set [T3] as input. The NLO QCD predictions 
have been taken from the plot in Fig. 3b) of [Tlj; these were obtained with the DISENT program 
[T2*] using MRST99 partons a renormalisation scale of Et and a factorisation scale of Q. 
The statistical, systematic and jet-energy- scale uncertainties have been added in quadrature 
to estimate the total experimental uncertainty. All the theory predictions give a reasonably 
good description of the data. The NLO predictions generally give a slightly better description 
than the LO predictions. For the (z, /^-factorisation approach, the "NLO" corrections are only 
significant in the forward region. 

In order to verify that the extra z convolution of (z, fc 4 )-factorisation with respect to k t - 
factorisation is important, we also repeated the calculation taking the limit z — * in the 
partonic cross section. (Usually, /c t -factorisation is only applied using the unintegrated gluon, 
whereas here we also include the unintegrated quark.) The parton emitted in the last evolution 
step then goes in the proton direction and is not counted in the inclusive jet cross section. In 
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125 < Q 2 < 10 5 GeV 2 and 8 < E T < 100 GeV 

i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 
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Figure 10: Comparison with ZEUS inclusive jet production data [H] at high Q 2 . The feint and 
bold lines correspond, respectively, to the predictions of the conventional QCD approach and 
the (z, /c t )-factorisation approach based on doubly-unintegrated parton distributions. 
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general, the predictions are much worse, even in the current jet region, providing evidence that 
the extra z convolution of our method is important. 



6.2 Comparison with HI data at low Q 2 

The HI Collaboration have measured the inclusive jet cross section in DIS at high Q 2 J7| and 
at low Q 2 ^H]- Here, we focus on the latter, where Q 2 = 5 to 100 GeV 2 . In this region, the NLO 
QCD corrections to LO QCD are larger than at high Q 2 , and the advantages of our approach 
become more apparent. Again, this data was taken during 1996 and 1997. 

The HI Collaboration use the electron method to reconstruct the kinematical variables, so 
cuts are imposed directly on the variable y, namely 0.2 < y < 0.6. We therefore set y min = 0.2 
and y m ax = 0.6 in ([50)1 . Also, HI present their data in rapidity bins in the lab frame rather 
than the Breit frame. It can be shown that the rapidity in the lab frame is 



vr = \ log 



AE ^ x f 2 \i + Xb ) ~ Xb ~ w 2et3tCOs4>ej ) f 



(98) 

where, in the Breit frame, the initial positron has transverse momentum squared of 



e 2 = AE p E e x B (AE p E e - 1) (99) 

and the 4-momentum of the outgoing jet has been written in the form 

j = cijp + bj q' + j±. (100) 

It is necessary to average the cross section over the azimuthal angle ej - between the positron 
and the outgoing jet in the transverse plane. For the "NLO" (z, fc^-factorisation calculation, 
the jet 4-momenta are not necessarily the same as the 4-momenta of the outgoing partons. It 
is necessary to pass the 4-momenta through a jet algorithm. Rather than use (|9*Hj) to determine 
77 LAB , which would require an additional azimuthal averaging, it is simpler to explicitly transform 
the 4-momenta from the Breit to the lab frame, then calculate the rapidity of the resultant 
4-momenta. 

In Fig. ^2 we show da/dE? integrated over Q 2 between 5 and 100 GeV 2 in three rapidity 
intervals. For the results presented, we used the MRST2001 LO parton set [THj as input. The 
NLO QCD predictions have been taken from the plot in Fig. 1 of the HI paper ^H]; these were 
obtained with the DISENT program using CTEQ5M partons JH], a renormalisation scale of 
Et and a factorisation scale of Q. The hadronisation correction factors used in jTHj have been 
applied to all the theory predictions. The statistical and systematic uncertainties have been 
added in quadrature to estimate the total experimental uncertainty. 

The LO (z, fcj)-factorisation calculation is in excellent agreement with conventional LO 
QCD, but neither describe the data well, especially in the forward rapidity region. The "NLO" 
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Figure 11: Comparison with HI inclusive jet production data [TH] at low Q 2 . The predictions of 
the (z, /ct)-factorisation approach based on doubly-unintegrated partons (which is much simpler 
to implement) are in good agreement with the conventional QCD approach. In some bins the 
predictions of the latter approach are hidden beneath the bold lines of the (z, /c t )-factorisation 
approach, at the respective order. 



29 



(z, /c t )-factorisation calculation is in very good agreement with conventional NLO QCD, al- 
though the agreement gets slightly worse as Et increases. 5 Deviations of the data from NLO 
QCD are seen only at small Et in the forward region. Here, the NLO corrections are quite large 
and it is likely that NNLO corrections or resolved virtual photon contributions are important 
in this region. Again, taking the limit z — > makes the (z, k t )-factorisation predictions much 
worse, showing that it is important to keep the precise kinematics. 

7 Conclusions 

We have presented a method for determining unintegrated parton distributions, f a (x, k 2 , (i 2 ), 
from the conventional (integrated) parton distributions, by considering the last DGLAP evolu- 
tion step separately, and imposing angular-ordering constraints on gluon emission. To include 
the precise kinematics in the hard subprocess initiated by the final parton in the evolution 
ladder, it is necessary to consider G?oufe%-unintegrated parton distributions, f a (x, z, k 2 , fi 2 ). We 
gave a prescription, called (z, /^-factorisation, for the computation of cross sections using these 
distributions. This prescription is a natural generalisation of the /^-factorisation approach. 

We used (z, /^-factorisation to estimate the cross section for inclusive jet production at 
HERA at lowest order. Using the same LO doubly-unintegrated distributions, we then carried 
out a "NLO" calculation which included the dominant Feynman diagrams with the soft gluon 
singularities being regulated by angular ordering. 

We showed that at O(ot s ) the predictions of the approach based on doubly-unintegrated 
partons, with exact kinematics, are close to the conventional LO QCD calculation at 0(a l s ). 
The relative simplicity of the former approach is shown schematically in Fig. ^ Similarly, at 
0{a\) the predictions of the approach based on doubly-unintegrated partons are close to the 
conventional NLO QCD calculation at 0{a 2 s ). 

It was seen that the NLO corrections are large in the forward region at low Et and low Q 2 
where the agreement with the data is poor. It is possible that the simplified (z, /c t )-factorisation 
approach might help to evaluate the role of the NNLO contribution. Alternatively, the resolved 
photon contribution is known to be important in the regime where Et is much greater than 
Q. It would be better to calculate the resolved photon contribution in terms of the doubly- 
unintegrated parton distributions of the photon. 6 

The logical next step would be to show that doubly-unintegrated distributions can be ap- 
plied to pp and pp collisions. The simplest calculation would be the transverse momentum 
distribution of produced W and Z bosons. 

5 In two bins the "NLO" (z, ^-factorisation predictions are significantly higher than the NLO QCD predic- 
tions. This is due to the jet algorithm applied, which increases the "NLO" (z, fct)-factorisation predictions by 
more than a factor of two in these two bins only, compared to the result when no jet algorithm is applied. 

6 In Ref. ^Hl) f° r example, the KMR prescription was applied to obtain the unintegrated gluon distribution 
of the photon. 
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We conclude that by reorganising the perturbative expansion in as to keep only the most 
important terms, our method provides a simple but effective way of estimating exclusive (and 
inclusive) observables. 
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